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Abstract 

> : 

^ ' The dynamics of the random-phase sine-Gordon model, which describes 2D 

vortex- glass arrays and crystalline surfaces on disordered substrates, is in- 
vestigated using the self-consistent Hartree approximation. The fluctuation- 
dissipation theorem is violated below the critical temperature T c for large time 
t > t* where t* diverges in the thermodynamic limit. While above T c the av- 
eraged autocorrelation function diverges as T\n(t), for T < T c it approaches 
a finite value q* ~ 1/(T C - T) as q(t) = q* - c{t/t*)~ v (for t -> i*) where z/ 
is a temperature-dependent exponent. On larger time scales t > t* the dy- 
namics becomes non-ergodic. The static correlations behave as ~ Tln|x| for 
T > T c and for T < T c when x < with £* ~ exp {A/(T C — T)}. For scales 

1 



x > £*, they behave as ~ m~ l T In \x\ where m T/T c near T c , in general 

agreement with the variational replica-symmetry breaking approach and with 

recent simulations of the disordered-substrate surface. For strong-coupling 

the transition becomes first-order. 
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Much attention has been given recently to the 2D sine-Gordon model []T] p~6|] with random 
phases which describes two very interesting disordered systems: (i) An array of flux-lines 
as obtained when a magnetic field is applied parallel to a type II superconducting film ||, 
and (ii) A 3D crystalline surface with a disordered substrate 0. It is also equivalent to a 
vortex-free XY model with a random field. 

The intensive theoretical attention is due to diverging results obtained by the renor- 
malization group (RG) [HHU, the Gaussian variational (GV) approximation n — > 

Bethe ansatz §|,|l(||, and by two extensive numerical studies |12|,[13|]. The Hamiltonian of the 



random sine-Gordon model (RSGN) is: 

H = J dx(^[Vh(x,t)} 2 - g cos (j[h(x,t) - . (1) 

The coefficient k is the stiffness, g is the coupling constant, 7 determines the periodicity 
(we choose 7 = 2n), and d(x) is a random variable such that jd(x) is a random phase 
uniformly distributed in (0, 2iv]. Random phases for different x are uncorrelated (short- 
distance correlations will not affect the conclusions). The variable h(x, t) denotes the vortex- 
line displacement from a periodic lattice or the height of the surface. 

The controversial issues concern the properties of the low-temperature phase. The quan- 
tity of interest is the correlation function: 

C{x -x',t- t') = (jh(S, t) - h{x', t')} 2 ) , (2) 

where (• • •) denotes averaging over the thermal noise and the quenched disorder. 

In the static limit, RG predicts that C(x — x') changes behavior from (T/ttk) In \x — 
x'\ for T > T c , to a behavior ~ e 2 (ln|x — x'\) 2 for T < T c where T c = k/tt is the critical 
temperature and e — \1 — T/T c \. The GV approach yields a one-step replica symme- 
try breaking (RSB) which gives correlations which remain logarithmic but with a different 
coefficient S(T) In \x — x'\ for T < T c . One study finds S(T) = T c /ttk while the other |7| 
finds a coefficient which increases from this value for temperature below T c . 

A potential resolution of the contradiction between the RG and the GV results was 



suggested recently by Le Doussal and Giamarchi . They have shown that at the replica- 



symmetric RG fixed-line for T < T c RSB is a relevant perturbation. This suggest the 
possibility of a spontaneous replica-symmetry breaking. 

Numerical simulations of the weak coupling RSGM show no sign of a transition what- 
soever in the equilibrium properties [|12|| . Our simulations of the random-substrate model, 
which is better described by the strong-coupling RSGM, exhibit a continuous transition at 
the analytically predicted temperature ||13|| . The behavior of the correlations is consistent 
with the predictions of the GV theory with one-step RSB. 

For the dynamical behavior the only analytic predictions are from heuristic arguments 
T5| , |T6H and from RG calculations |J. The RG predicts anomalous dynamics for T < T c which 



is manifested by a temperature-dependent dynamic exponent such that z = 2(1 + 1.78e) and 
a non-linear response to a small driving force F. The average velocity is predicted to be 
proportional to F z l 2 for small F. The simulations |]12| vindicate qualitatively this latter 



prediction although the coefficient of e is much smaller than that predicted by dynamic RG. 

In view of these results we study here a non-perturbative self-consistent approach to 
the dynamics of the system as a complement to the RG calculations. The general ap- 
proach we take is based on that introduced by Sompolinsky and Zippelius (SZ) |T7| for the 
Sherrington-Kirkpatrick (SK) model. This approach has been extended by Crisanti, Horner, 
and Sommers for the spherical p-spin interacting spin-glass model |18| and generalized by 
Kinzelbach and Horner (KH) [19| for fluctuating manifolds in disordered media. 

We begin by the RSGM Hamiltonian ([!]) with random-quenched phase with zero mean 
and correlations: 

e Vid Me- vld & ) ) = 5{$-£). (3) 
The stochastic equation of motion obeyed by h(x, t) within purely relaxational dynamics is: 

1 »*<*.*>- . p "m_ + (i M. (4 ) 



T dt 8h(x, t) 

The Gaussian white noise C(^>*) obeys (((x,t)((x',t')) = (2/T )<5(x — x')5(t — t'). T is a 
microscopic kinetic coefficient and (3 = 1/T. The next step is to obtain in a standard way 



the Martin-Siggia-Rose (MSR) generating functional |^TJ which may directly be averaged 
over the noise ( and the quenched disorder d(x). The action so obtained reads: 



S = I dx I dt\-Yl 



ih(x, t) + ih(x, t) d t - T k(3V 2 h(x, t) 



- -f3 2 T 2 ^ 2 g 2 I dx I dtdt' ih(x, t)ih(x, t') cos [7 (h(x, t) - h(x, £'))] • 



(5) 



In the self-consistent Hartree approximation all terms which contains more than two 
fields h or h are written in all possible ways as products of such pairs and all pairs but 
one are replaced by their expectation values. The action then is quadratic in h and h with 
coefficients which are only functions of response and correlation functions, defined as: 



r(x - f , t - t') = (h(x, t)zh(x', £')) , t>t' 
q(x - f , t - t') = (h(x, t)h(x, t)) - (h(x, t)h(x', t')) . 



(6) 
(7) 



Their Fourier coefficients are denoted by r(k,t) and q(k,t). The effective action obtained in 
this way is: 



S e . 



dx dt< -r 



ih(x,t) +ih(x,t) d t - T KpV 2 + fi(0) h(x,t) 



1 



dx I dtdt' \ -u(t - t')ih(x, t) ih(x, t') + X(t - t')ih(x, t)h(x, f) \ . (8) 



where we have defined (r = 1,7 = 2-n) : 
1 



u(t) = ^ 2 Tfa 2 g 2 (cos [2tt (h(x , t) - h(x, 0)]> r . 



(9) 



X(t) = 4:7i 2 uj(t)r(t) and /i{t) = ff°dt'\(t'). The functions lu, X and fi depend on t only 
through q(t) = q(0,t) and r(t) = r(0,t). 

The next step is to obtain the equations of motion. By standard techniques explained 
in details in Ref. |T9|], in momentum space they read: 



{d t + (3Kk 2 + n(o)}r(k,t) - j dt'X{t - t')r(k,t') =0 



(10) 



with initial condition r(fc,0 + ) = 1, r(k,t) = for t < 0, and: 
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[d t + (3nk 2 + fj,(o)}q{k,t) - ! dt'X(t-t')q(k,t') = 1 - {l(k, t) - I(k, 0)} (11) 



Jo 



with q(k,0) = and I(k,t) = / °° dt' {u(t + f)r(k, t') - X(t + t')q(k, t')\. If fluctuation- 
dissipation theorem (FDT) holds, the I(k, t) — I(k, 0) = and equations for r and q become 



equivalent as they should. The validity of the FDT implies that for long time: 



(3nk 2 + n{t)] q(k,t) < 1. 



(12) 



This condition may be cast as: 



?(*)< 



Ann 



T 



In 1 + 



n(t)T 



) 



(13) 



Following KH [|T^] we define the function: 



A(?) 



i 2g 2 




(14) 



T e ^ K il T - 1 kT 



in terms of which the FDT dynamics holds as long as A(g) > 0. However, this condition 
will be violated at low enough temperatures. For g > g tr = k/\/2tx the transition is "first 
order" in the sense that as T is lowered A(g) vanishes at a finite q. We shall postpone 
the discussion of the first-order regime and the tricritical behavior which takes place at gt r 
to a future, longer, publication. Here we concentrate on the regime g < g tr for which the 
transition is continuous. In this regime A(q) becomes negative at the (7-independent critical 
temperature T c = k/tt (as found by other analytic approaches). In this regime, and close to 
T c , A(q) < for q > q x = (l/47r 2 e)| ln(2# 2 / ' kT)\. 

The solution proposed assumes the existance of a sharp separation of time scales between 
the ergodic behavior (i.e., relaxation within pure states) for t <C t* and the non-ergodic 
behavior for t 3> t* which includes "tunneling" between different pure-states. In this latter 
regime X(t) < \d t u\ and the integral fj?dt'X(t') = [i* > 0. Therefore for t < £*, A(q) = 
A(q)+/j l * and the instability is avoided. Within the so-called "quasi-FDT" (QFDT) approach 
|T8|420[ it is assumed that \x* is the minimal one necessary to make A non-negative. This 
choice implies A'(q*) = 0, where q* is derived from the ergodic stationary solution as t — > t*\ 
q*(k) = l/((3Kk 2 + //*) and therefore q* = (T/4ttk) ln(l + ktt 2 /h*T). 
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At long-time t >f the time variable rescales to r = t/t*. We look for the solutions for 
the response and correlation functions which have the scaling forms q(k,t) = q(k,t/t*) and 
r(k,t) = (l/t*)r(k,t/t*) with similar expressions for u>(t) and A(t). In the limit t —>■ oo the 
equations for q and f have solution of the form: 



md T q(k, r) = r(k, r) 



t > t* 



(15) 



which implies: 



md T uj{r) = — A(t), 



t > t*. 



(16) 



Then the equations for q and r are equivalent and hence quasi-FDT holds. Equation (|16"D 
also determines m since its integration leads to: 



m 



(17) 



where u* = u)(q*). Using the equation for q*, (A'(q*) = 0), we find 

TFT 



m(T) 



1 _ e -4™?7T 



(18) 



We note that as t — > t* the solution for t > t* has to merge with the one below t*. This 
matching condition yields the equation: 

rffc 1 



1 = 4ttV 



(19) 



(27r) 2 (/3 K A; 2 + y u*) 2 

which is exactly the condition A'(q*) = 0, or in the replica language it is the condition on 
/i* for the eigenvalue of the replicon mode to change from a negative value to zero (see Eq. 
(5) in Ref. §). 

For t ^> t* — > oo, the stationary value of g(/c) is now given by: 

T 



qoo(k) 



nk 2 



which leads to: 



<?oo(fc) 



rank 2 



1 + (1 -m)u)*q*(k) 



1 — m T 

m K,(k 2 + h*T/k) 



(20) 



(21) 



In the regime T <T C the equal-time correlations are: 

— In \x\ if x <C £ * 

C{x) = \ «" (22) 

(l - e - 4wKq */ T ) In |f | if x > £* 



7TK 



where £* = y k/T/i* ~ exp(A/e). As T c is approached from below the slope of the logarithm 
of the asymptotic behavior is given by: 

1 



7l 2 



1 + 



2V/k 2 1/6 ) (23) 



which is consistent with the behavior found by Korshunov 0. Fig. [I] shows the slope S(T) 
for different values of g. Near T c the behavior is also consistent with Ref. M and with the 



numerical results we obtained simulating the random-substrate model [13 . 

The dynamics within each ergodic component for time t < t* for T < T c is studied by 
assuming a solution of the form 

q{t) = q*-6q{t). (24) 
Expanding in 8q(t) we find a solution of the form 8q(t) ~ t~ u where v is a temperature 



dependent exponent determined from the solution of the equation [jTB[|19 

Y 2 (l-u 



m(T), (25) 



T(l-2u) 

where T is the gamma function. The plot of v{T) is given in Fig. ||. 

So, while above T c , q(t) diverges as (T/47m) ln(t), below T c it approaches a finite value 
as t — > t* (and t* — > oo as L — > oo ), g(t) = g* - c As T — > T~, z/ — >■ and 



the approximation (p4 ) should be handled carefully. The asymptotic behavior becomes 
logarithmically slow. In this limit the value of q* diverges like q* ~ 1/ (T c — T) . 

If we keep the size L dependence explicitly and take the limit L — > oo after t — > oo, the 
behavior of q(t) will be modified as T c is approached from above. This will be studied in a 
future paper, where we shall also address the first order transition which occurs in the strong 
coupling regime. In this regime the dynamic critical temperature is g-dependent and larger 



S 



than the static one. Such a behavior was found in other systems with first-order transition 
I8| , [19|j and it signals the existence of many metastable states ||22j| . 

To summarize, in the second-order regime the self-consistent dynamic approach yields 
the following results: In equlibrium it essentially reproduces the results of the GV approach 
with the one-step replica symmetry breaking. In the dynamics this is refelected by the 
breaking of the FDT below T c on large time scales t > t*. Phase space is broken into 
separate components and the divergence of their barriers leads to the divergence of t* in the 
thermodynamic limit. Within each component the autocorrelations q(t) approach a finite 
value q*. At short time q(t) ~ hi(t) and q(t) ~ q* — ct~ u as t —>■ t*. q* diverges and v — ► 
as T — > T~. This behavior is obtained based on the self-consistent quasi-FDT which holds 
for t > t*. The analysis applies if t — > oo before L — > oo, where L is the system size. If the 
order of limits is reversed the system never reaches equilibrium below T c . This will entail a 
non-stationary behavior with aging phenomena |23| . 

This research was supported in part by the National Science Foundation under Grant 
No. PHY89-04035. Acknowledgment is also made to the donors of The Petroleum Research 
Fund, administrated by ACS, for partial support of this research. 
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FIGURES 

FIG. 1. The slope of the ln|x| term in Eq. (^) versus temperature for different coupling 
constants g < gt r and k = 1. 

FIG. 2. The temperature dependence of the critical exponent v for different coupling constants 
g. As g decreases they approach the solution of the Eq. (|24|) with m = ttT/k. 
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